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ABSTRACT 

We report the discovery of HAT-P-24b, a transiting extrasolar planet orbiting the moderately bright 
V=11.818 F8 dwarf star GSC 0774-01441, with a period P = 3.3552464 ± 0.0000071 d, transit epoch 
T c = 2455216.97669 ± 0.00024 (BJE0), and transit duration 3.653 ± 0.025 hours. The host star has a 
mass of 1.191 ± 0.042M o , radius of 1.317 ± 0.068 R Q , effective temperature 6373 ± 80K, and a low 
metallicity of [Fe/H] — —0.16 ± 0.08. The planetary companion has a mass of 0.681 ± 0.031 Mj, and 
radius of 1.243 ± 0.072 i?j yielding a mean density of 0.439 ± 0.069gcm~ 3 . By repeating our global 
fits with different parameter sets, we have performed a critical investigation of the fitting techniques 
used for previous HAT planetary discoveries. We find that the system properties are robust against 
the choice of priors. The effects of fixed versus fitted limb darkening are also examined. HAT-P-24b 
probably maintains a small eccentricity of e = 0.0521q'q^, which is accepted over the circular orbit 
model with false alarm probability 5.8%. In the absence of eccentricity pumping, this result suggests 
HAT-P-24b experiences less tidal dissipation than Jupiter. Due to relatively rapid stellar rotation, 
we estimate that HAT-P-24b should exhibit one of the largest known Rossiter-McLaughlin effect 
amplitudes for an exoplanet (AVrm — 95m/s) and thus a precise measurement of the sky-projected 
spin-orbit alignment should be possible. 

Subject headings: planetary systems — stars: individual (HAT-P-24, GSC 0774-01441) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

The understanding of planetary systems has been 
spearheaded by the study of transiting extrasolar planets 
in recent years. Eclipses have long offered a key to un- 
locking the secrets of the heavenly bodies, for example in 
Solar System studies and in the field of eclipsing binaries. 
For exoplanets, an eclipse offers a door into the inner- 
workings of an alien system hundreds of light-years away. 
As an inherently low probability event, each and every 
transiting system is precious and timeless to the plane- 
tary scientist. Transits continue to offer unprecedented 
access to an exoplanet's nature and allow for, amongst 
other things, the determination of the oblateness of a 
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planet (|Seaeer fc Huil (|2002T) : ICarter fc Wirml mW\). 
therm al mapping of the planetary surface (Knutson et al. 
2007) and accurate plane tary radii at the percent level 
( Charbon neau et al.N2000D . 

The Hu ngarian-made Aut omated Telescope Network 
(HATNet: IBakos et aT1 l2004) survey has been one of the 
principal contributors to the discovery of transiting ex- 
trasolar planets (TEPs). In operation since 2003, it has 
now covered approximately 14% of the sky, searching 
for TEPs around bright stars (8 < / < 14.0). HAT- 
Net operates six wide-field instruments: four at the Fred 
Lawrence Whipple Observatory (FLWO) in Arizona, and 
two on the roof of the hangar housing the Smithso- 
nian Astrophysical Observatory's Submillimeter Array, 
in Hawaii. Since 2006, HATNet has discovered 23 TEPs 
(with 16 announced or published so far). In this work we 
report our 24th discovery, around the relatively bright 
star previously known as GSC 0774-01441. 

In Section [2]we report the detection of the photometric 
signal and the follow-up spectroscopic and photometric 
observations of HAT-P-24. In Section [3] we describe the 
analysis of the data, beginning with the determination 
of the stellar parameters, continuing with a discussion of 
the methods used to rule out non-planetary, false positive 
scenarios which could mimic the photometric and spec- 
troscopic observations, and finishing with a description 
of our global modeling of the photometry and radial ve- 
locities. In Section 2] we investigate the variations of our 
results using different parameter sets with uniform priors 
to test the robustness of the fitted parameters. We also 
discuss the effects of fitting versus fixing limb darkening 
coefficients. In Section [5] to Section [SJ we present discus- 
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sions and analyses of the orbital eccentricity, a possible 
linear drift in the RVs, the measured mid-transit times 
and possibilities for future follow-up observations. Fi- 
nally, we summarize our findings in Section [9] 

2. OBSERVATIONS 

2.1. Photometric detection 

The transits of HAT-P-24b were detected with the 
HAT-5 and HAT-6 telescopes in Arizona, and with the 
HAT-8 and HAT-9 telescopes in Hawaii. The star 
GSC 0774-01441 lies in the overlap of two fields, in- 
ternally labeled as 315 (07:30 +15:00) and 314 (07:00 
+15:00). The former field was observed on a nightly ba- 
sis between 2007 October and 2008 May, while the lat- 
ter field was observed between 2008 November and 2009 
May. For field 315 we gathered 8551 R band exposures of 
5 minutes at a 5.5 minute cadence, while for field 314 we 
gathered 5503 Sloan r band images with the same expo- 
sure time and cadence. Each field 315 image contained 
approximately 51,000 stars down to R ~ 14, while each 
314 image contained approximately 130,000 stars down 
to r ~ 14.5. For the brightest stars in field 315, we 
achieved a per-image photometric precision of 3 mmag, 
while for field 314 we achieved 5 mmag precision. 
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Fig. 1.— Unbinned light curve of HAT-P-24 including all 14,000 
instrumental R band and Sloan r band 5.5 minute cadence mea- 
surements obtained with the HAT-5, HAT-6, HAT-8 and HAT-9 
telescopes of HATNet (see the text for details), and folded with the 
period P = 3.3552401 days resulting from the global fit described 
in Section[3|l. The solid li ne sh ows the "P1P3" transit model fit 
to the light curve (Section 13.31 1. Solid squares show the 20-point 
binned light curve. 

The calibration of the HATNet frames was carried out 
using standard photometric procedures. For field 314 
the calibrated images were then su bjected to star detec - 
tion and astrometry, as described in lPal fc Bakosl (|2006(1 . 
Aperture photometry was performed on each image at 
the stellar centroids derived from the Tw o Micron All 
Sky Survey (2MASS: ISkrutskie et ail 12006( 1 catalog and 
the individual astrometric solutions. For field 315 we 
performed image su btraction ph otometry following the 
methods described in[Pal (2009b). Aperture photometry 
was performed on the difference images at the stellar cen- 
troid derived from 2MASS, while for the reference flux 
we adopted the r magnitude of each star, transformed 
from its 2MASS J, H and K$ magnitudes, and made 
use of the average relation between r and flux measured 
on the reference image via aperture photometry. The 
resulting light curves for both fields were decorrelated 



(cleaned of trends) using the Externa l Parameter Decor- 
relation (EPD; see B akos et al.l 12010( 1 technique in "con- 
sta nt" mode and the Trend Filtering Algorithm (TEA; 
see IKovacs et a\W 005). The light curves from each field 
were independently searched for periodic b ox-shaped sig- 
nals u sing the Box Least-Squares (BLS; see IKovacs et al.1 
2002) method. We detected a significant signal in both of 
the HATNet light curves of GSC 0774-01441 (also known 
as 2MASS 07151801+1415453; a = 07 h 15 m 18.00s, 5 = 
+14°15'45.4"; J2000; V=11.818 IDroege et all l200fih . 
with an apparent depth of ~ 7.2 mmag, and a period 
of P = 3.3552 days (see Figured]). 

2.2. Reconnaissance Spectroscopy 

As is routine in the HATNet project, all candidates are 
observed spectroscopically initially to establish whether 
the transit-like feature in the light curve is of non- 
planetary origin such a grazing eclipsing binary (i.e. a 
false positive). For example, large radial- velocity varia- 
tions of the star (tens of kms -1 ) would indicate such a 
circumstance. 

To perform this task, we used the Harvard- 
Smithsonian Cent er for Astrop hysics (CfA) Digital 
Speedometer (DS; Latham] I1992H : an echelle spectro- 
graph mounted on the FLWO 1.5 m telescope. This 
instrument delivers high- resolution spectra (A/ A A w 
35,000) over a single order centered on the Mg lb triplet 
(~5187A), with typically low signal-to-noise (S/N~ 10) 
ratios that are nevertheless sufficient to derive radial ve- 
locities (RVs) with moderate precisions of 0.5-1. Okms -1 
for slowly rotating stars. The same spectra can be used 
to estimate the effective temperature, surface gravity, 
and projec ted rotationa l veloc ity of the host star, as de- 
scribed bv lTorres etail (|2002f l. With this facility we are 
able to reject many types of false positives, such as F 
dwarfs orbited by M dwarfs, grazing eclipsing binaries, 
or triple or quadruple star systems. 

For HAT-P-24 we obtained five observations with the 
DS between 2008 May and 2009 January. The velocity 
measurements showed an r.m.s. residual of 0.74km.s , 
consistent with no detectable RV variation within the 
precision of the measurements. All spectra were single- 
lined, i.e., there is no evidence for additional stars in 
the system. The atmospheric parameters we infer from 
these observations are the following: effective temper- 
ature T g* = 7000 ± 100 if, surface gravity logg* = 
4.5 ± 0.25 (log cgs), and projected rotational velocity 
vsmi = 11.2 ± l.Okms -1 . The effective temperature 
corresponds to a F8 dwarf. The mean heliocentric RV of 
HAT-P-24 is 7rv = -2.09 ± 0.33 kms" 1 . We stress that 
the DS stellar parameters are based on solar composition 
models. 

2.3. High resolution, high S/N spectroscopy 

Given the significant transit detection by HATNet, 
and the encouraging DS results that rule out obvious 
false positives, we proceeded with the follow-up of this 
candidate by obtaining high-resolution, high-S/N spec- 
tra to characterize the RV variations, and to refine the 
determination of the stellar parameters. For this we 
used the HIRES instrument ()Vogt et al.l I1994H on the 
Keck I telescope located on Mauna Kea, Hawaii, between 
2009 April and 2009 December. The width of the spec- 
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trometer slit was 0'.'86, resulting in a resolving power of 
A/AA « 55,000, with a wavelength coverage of ^3800- 
8000 A. 

We obtained 18 exposures through an iodine gas ab- 
sorption cell, which was used to superimpose a dense 
forest of I2 lines on the stellar sp ectrum and establish a n 
accurate wavelength fiducial (see lMarcv fc Butler]|1992t) . 
An additional two exposures were taken without the io- 
dine cell, for use as templates in the reductions. In 
practice, we used only the second higher S/N template 
spectrum. Relative RVs in the sola r system barycentri c 
frame were derived as described bv lButler et al.l (|1996l ). 
incorporating full modeling of the spatial and temporal 
variations of the instrumental profile. The RV measure- 
ments and their uncertainties arc listed in Table [T] The 
period-folded data, along with a best fit described below 
in Section El are displayed in Figure [2j 

150 | 1 1 1 , 1 1 1 
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matching each spectrum to a reference spectrum using 
a transformation that includes a wavelength shift and a 
flux scaling that is a polynomial as a function of wave- 
length. The transformation was determined on regions 
of the spectra that are not used in computing this in- 
dicator. Note that our rel ative S index has not been 
calibr ated to the scale of IVaughan. Preston fc Wilsonl 
(|1978ft . The RMS of the relative S values is 3%, which 
is higher than the median formal error of 0.6% based 
on photon statistics, however the errors in this case are 
likely dominated by systematics in the spectrum match- 
ing procedure, which are difficult to quantify, so we do 
not consider this to be a robust detection of variability. 
We note that a 3% variati on is comparable to t hat found 
for other late F stars (e.g. iShkolnik et al.ll20"08l measured 
a ~ 1% median absolute deviation in the fluxes of the K 
line cores for r Boo and HD 179949, which corresponds 
to a similar expected RMS). 

We also note that S is uncorrelated with orbital phase; 
such a correlation might have indicated that the RV vari- 
ations could be due to stellar activity, casting doubt on 
the planetary nature of the candidate. There is no sign 
of emission in the cores of the Ca II H and K lines 
(S/N~ 38) in any of our spectra, from which we con- 
clude that the chromospheric activity level in HAT-P-24 
is very low. 

2.4. Photometric follow-up observations 
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Fig. 2. — Top panel: Keck/HIRES RV measurements for 
HAT-P-24 shown as a function of orbital phase, along with our 
best-fit eccentric orbit model (solid) and circular orbit model 
(dashed). Zero phase corresponds to the time of mid-transit. The 
center-of-mass velocity has been subtracted. Second panel: Phase 
folded velocity O — C residuals from the best fit. The error bars 
include a component from astrophysical jitter (7.4 ms" 1 ) added 
in quadrature to the formal errors (see Section l3.3jl . Third panel: 
Bisector spans (BS), with the mean value subtracted. The mea- 
surement from the template spectrum is included (see Section 13.211 . 
Bottom panel: Relative chromospheric activity index S measured 
from the Keck spectra. Note the different vertical scales of the 
panels. 

In the same figure we show also the relative S index, 
which is a measure of the chromospheric activity of the 
star derived from the flux in the cores of the Ca II H and 
K lines. This i ndex was computed following the pr escrip- 
tion given by IVaughan. Preston fc Wilsonl (|1978l ) , after 
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Fig. 3. — Unbinned normalized Sloan i band transit light curves, 
acquired with KeplerCam at the FLWO 1.2 m telescope on 2010 
Jan 09, 2010 Jan 16, and 2010 Feb 05, and with the Faulkes Tele- 
scope North on 2010 Feb 15. The li ght c urves have been EPD 
and TFA processed, as described in § 13.31 Final fits come from 
the {p 2 ,6,Ti.5 3.5} light curve parameter set and fitted limb dark- 
ening, in combination with an eccentric orbit RV fit. The dates of 
the events are indicated. Curves after the first are displaced ver- 
tically for clar ity. Our best fit from the global modeling described 
in Section 13.31 is shown by the solid lines. Residuals from the fits 
are displayed at the bottom, in the same order as the top curves. 
The error bars represent the photon and background shot noise, 
plus the readout noise. 

In order to permit a more accurate modeling of the 
light curve, we conducted additional photometric ob- 
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TABLE 1 

Relative radial velocities, bisector spans, and activity index 
measurements of hat-p-24. 



BJD 


RV a 


o-RV b 


BS 


CBS 


S c 


<^s c 


(2,454,000+) 


(ms^ 1 ) 


(m s _1 ) 


(ms- 1 ) 


(ms- 1 ) 






928.76217 . 


—49.49 


6.43 


1.10 


1.03 


0.44 


0.003 


954.79832 . . 






—4.52 


1.35 


0.43 


0.003 


9oo.oU587 . . 


—85.15 


7.37 


—47.96 


4.39 


0.42 


0.008 


flCf? Q1 ,4 fl ,4 

yoo.ol4y4 . . 


— 7.84 


6.61 


—65.72 


5.91 


0.45 


0.005 


nco with 

963.78272 . . 


50.62 


6.41 


2.42 


0.68 


0.45 


0.003 


1 1 rtf7 1 /I COO 

1107.14583 . 


—80.20 


8.64 


25.58 


1.21 


0.45 


0.003 


1109.0656 f . 


48.40 


7.00 


8.95 


0.13 


0.44 


0.003 


1112.13808 . 


84.68 


9.02 


12.27 


0.08 


0.43 


0.003 


1134.10472 . 


-78.19 


6.49 


1.68 


0.92 


0.45 


0.003 


1135.08419 . 


55.31 


5.83 


8.32 


0.41 


0.44 


0.002 


1136.10198 . 


65.26 


7.53 


1.35 


0.77 


0.43 


0.002 


1172.94629 . 


51.15 


6.80 


16.94 


0.21 


0.43 


0.002 


1174.10042 . 


-85.04 


6.32 


7.30 


0.48 


0.43 


0.001 


1187.96722 . 


-79.11 


7.17 


7.57 


0.35 


0.40 


0.002 


1189.07232 . 


71.01 


5.53 


1.20 


0.80 


0.43 


0.001 


1191.03182 . 


-80.15 


5.68 


5.56 


0.56 


0.44 


0.002 


1192.06936 . 






2.08 


0.77 


0.44 


0.001 


1192.07902 . 


37.54 


5.33 


3.51 


0.66 


0.44 


0.001 


1193.04199 . 


49.79 


6.10 


5.20 


0.59 


0.43 


0.002 


1193.85228 . 


-60.65 


5.88 


7.09 


0.52 


0.43 


0.002 



Note. — Note that for the iodine-free template exposures we do not 
measure the RV but do measure the BS and S index. Such template 
exposures can be distinguished by the missing RV value. 
a The zero-point of these velociti es is arbitrary. An overall offset 7 re i 
fitted to these velocities in Section l3.3l has not been subtracted. 
b Internal errors excl uding the component of astrophysical jitter con- 
sidered in Section 13.31 

c Relative chromospheric activity index, not calibrated to the scale of 
IVaughan. Preston fc Wilson! lfl978l 1. 



servations with the KeplerCam CCD camera on the 
FLWO 1.2 m telescope in Arizona and with the 2.0 m 
Faulkes Telescope North (FTN) at Haleakala Observa- 
tory in Hawaii. We observed three transit events of HAT- 
P-24 with the FLWO 1.2 m telescope on the nights of 
2010 Jan 09, 2010 Jan 16 and 2010 Feb 05, and a fourth 
transit event with the FTN on the night of 2010 Feb 15 
(Figure [3J . These observations are summarized in Ta- 
ble n 

The reduction of these images, including basic cali- 
bration, astrometry, and aperture photomet ry, was per- 
formed as described by iBakos et al.l (|2010T ) . We per- 
formed EPD and TFA to remove trends simultaneously 
with the light curve modeling (for more details, see 
Section [3j and lBakos et al.ll2010D . The final time series 
are shown in the top portion of Figure G2 along with our 
best-fit transit light curve model described below; the 
individual measurements are reported in Table El 

3. ANALYSIS 

3.1. Properties of the parent star 

Fundamental parameters of the host star HAT-P-24 
such as the mass (M*) and radius (R*), which are 
needed to infer the planetary properties, depend strongly 
on other stellar quantities that can be derived spec- 
troscopically. For this we have relied on our template 
spectrum obtained with the Keck/HIRES instrument, 
and the anal ysis package known as Sp ectroscopy Made 
Easy (SME; IValenti &: Piskunovl I1996D. a l ong w ith the 
atomic line database of IValenti fc Fischerl (|2005l ). SME 
yielded the following initial values and uncertainties 
(which we have conservatively increased for T e ff* and 



[Fe/H] to include our estimates of the systematic er- 
rors): effective temperature T c ff* = 6188 ± 80 K, stel- 
lar surface gravity log g* = 4.01 ± 0.06 (cgs), metallicity 
[Fe/H] = —0.26 ± 0.08 dex, and projected rotational ve- 
locity i>sini = 10.5 ± 0.5kms _1 . 

In principle the effective temperature and metallicity, 
along with the surface gravity taken as a luminosity in- 
dicator, could be used as constraints to infer the stellar 
mass and radius by comparison with stellar evolution 
models. For planetary transits a stronger constraint is 
often provided by the a/i?* normalized semi- major axis, 
which is closely related to p+, the mean stellar density. 
The quantity a I R* can be derived directly from the tran- 
sit light curves (jSeager fc Mallen-Ornelas1l2003D and the 
RV data (for eccentric cases, see Kipping (|2010T l). This, 
in turn, allows us to improve on the determination of the 
spectroscopic parameters by supplying an indirect con- 
straint on the weakly determined spectroscopic value of 
log that removes d egeneracies. We tak e this approach 
here, as described in IBakos et a l. (2010). The validity 
of our assumption, namely that the adequate physical 
model describing our data is a planetary transit (as op- 
posed to a blend), is shown later in Section [3721 

After the first iteration for determi n ing t he stellar 
properties, as described in IBakos et al.l (|2010l ). we find 
that the surface gravity, logg* = 4.27 ± 0.04, is signifi- 
cantly different from our initial SME analysis, which is 
not surprising in view of the strong correlations among 
T G ff*, [Fe/H], and logg* that are often present in spec- 
troscopic determinations. Therefore, we carried out 
a second iteration in which we adopted this value of 
log<?* and held it fixed in a new SME analysis (cou- 
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TABLE 2 

Summary of photometric follow-up observations, all of which were 
taken in Sloan i band. 



Facility Date Number of Images Cadence (s) 

KeplerCam/FLWO 1.2 m 2010 Jan 09 527 39 

KeplerCam/FLWO 1.2 m 2010 Jan 16 246 53 

KeplerCam/FLWO 1.2 m 2010 Feb 05 256 39 

FTN 2010 Feb 15 534 30 



TABLE 3 

High-precision differential photometry of 
HAT-P-24 



BJD 


Mag a 


CMag Filter 


(2,400,000+) 






55206.63289 


-0.00145 


0.00176 i 


55206.63724 


0.00220 


0.00181 i 


55206.63790 


-0.00278 


0.00173 i 


55206.64440 


-0.00213 


0.00174 i 


55206.64485 


0.00093 


0.00171 i 


55206.64548 


-0.00051 


0.00161 i 


55206.64594 


0.00535 


0.00161 i 


55206.64659 


-0.00391 


0.00141 i 


55206.64703 


0.00162 


0.00140 i 


55206.64767 


-0.00023 


0.00134 i 



Note. — The complete table is available in a 
machine-readable form in the on-line journal. A por- 
tion is shown here for guidance regarding its form and 
content. 

a The out-of-transit level has been subtracted. These 
magnitudes have been subjected to the EPD and TFA 
procedures, carried out simultaneously with the transit 
fit. 

pled with a new global modeling of the RV and light 
curves), adjusting only T e g*, [Fe/H], and usini. This 
gave T cSir = 6373 ± 80 K, [Fe/H] = -0.16 ± 0.08, and 
wsini = 10.0 ± 0.5 kms , in which the conservative un- 
certainties for the first two have been increased by a fac- 
tor of two over their formal values, as before. Experience 
with the SME analysis for previous HAT planets leads us 
to select this estimate. A further iteration did not change 
log significantly, so we adopted the values stated above 
as the final atmospheric properties of the star. They are 
collected in Table HI together with the adopted values for 
the macroturbulent and microturbulent velocities. 

The low metallicity of HAT-P-24 is interesting in that 
a well-known bias exists for finding gia nt planets around 
metal- rich stars () Johnson et al.l 120101) . Uohnson et al.l 
(2010) find that the occurrence of giant planets scales 
as / ~ ]_Qi-2[FiB/H] anc j a j g0 re p 0r t a scaling with stellar 
mass of / ~ . This means that the a-priori probabil- 
ity of finding a planet around HAT-P-24 is a respectable 
~ 75% that of a Solar-like star. 

With the adopted spectroscopic parameters the model 
isochrones yield the stellar mass and radius M± = 1.191± 
0.042 Mq and R* = 1.317 ± 0.068 i?©, along with other 
properties listed at the bottom of Table |U HAT-P-24 is 
a F8 dwarf star with an es timated age of 2.8 ± 0.6 Gyr, 
according to these models (|Yi et al.l l2001). The inferred 
location of the star in a diagram of a/R* versus T e ff*, 
analogous to the classical H-R diagram, is shown in Fig- 
ure |U The stellar properties and their ler and 2<r con- 
fidence e l lipsoi ds are displayed against the backdrop of 
lYi et al.l (|2001l ) isochrones for the measured metallicity 
of [Fe/H] = —0.16, and a range of ages. For comparison, 



the location implied by the initial SME results is also 
shown (triangle), and corresponds to a somewhat more 
evolved state. 
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Fig. 4. — Model isochrones from Yi et al. (2001) for the measured 
metallicity of HAT-P-24, [Fe/H]= -0.16, and ages of 0.2, 0.5, 1.0, 
2.0, 3.0, 4.0, 5.0, 6.0, 7.0 and 8.0 Gyr (left to right). The adopted 
values of T e fj* and a/R* are shown together with their lcr and 
2cr confidence ellipsoids. The results come from the {p a ,6,Ti.5 s.g} 
light curve parameter set and fitted limb darkening, in combination 
with an eccentric orbit RV fit. The initial values of T e g* and a/H* 
from the first SME and light curve analyses are represented with 
a triangle. 

The stellar evolution modeling provides color indices 
that may be compared against the measured values as 
a consistency check. The best available measurements 
are the near-infrared magnit udes from the 2MASS Cat- 
alogue (|Skrutskie et al.ll2006h . J 2 mass = 10.797 ± 0.022, 
#2MASS = 10.589 ±0.024 and i^MASS = 10.543 ± 0.020; 
which we have converted to the photometric system of 
the models (ESQ system) using the transformations by 
(|Carpenterll2001l ). The resulting measured color index is 
J-K = 0.274± 0.032. This is within la of the predicted 
value from the isochrones of J—K = 0.29±0.02. The dis- 
tance to the object may be computed from the absolute 
K magnitude from the models (Mr = 2.59 ± 0.11) and 
the 2MASS K s magnitude, which has the advantage of 
being less affected by extinction than optical magnitudes. 
The result is 396 ± 20 pc, where the uncertainty excludes 
possible systematics in the model isochrones that are dif- 
ficult to quantify. 

3.2. Spectral line-bisector analysis 

Our initial spectroscopic analyses discussed in 
Section |2"T21 and Section |2~31 rule out the most obvious as- 
trophysical false positive scenarios. However, more sub- 
tle phenomena such as blends (contamination by an un- 
resolved eclipsing binary, whether in the background or 
associated with the target) can still mimic both the pho- 
tometric and spectroscopic signatures we see. 
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TABLE 4 
Stellar parameters for HAT-P-24 



Parameter Value Source 



Spectroscopic properties 

T eff * (K) 6373 ± 80 SME a 

[Fe/H] -0.16 ±0.08 SME 

vsini (kms" 1 ) . . . 10.0 ± 0.5 SME 

Umac (kms" 1 ).... 4.94 SME 

v mic (kms" 1 ) .... 0.85 SME 

7rv (kms" 1 ) -2.09±0.74 DS 

Photometric properties 

V (mag) 11.818 TASS 

V-I c (mag) 0.628 ±0.089 TASS 

J (mag) 10.797 ± 0.022 2MASS 

H (mag) 10.589 ± 0.024 2MASS 

K s (mag) 10.543 ± 0.020 2MASS 

Derived properties 

M 4 (Mq) 1.191 ±0.042 YY+a/i?*+SME b 

Ri, (R©) 1.317 ±0.068 YY+a/R*+SME 

logg* (cgs) 4.27 ±0.04 YY+a/iJ*+SME 

L* (Lq) 2.56 ±0.31 YY+a/i?*+SME 

M v (mag) 3.74 ±0.14 YY+a/.R*+SME 

M K (mag,ESO) . . 2.59 ±0.11 YY+a/i?*+SME 

Age (Gyr) 2.8 ± 0.6 YY+a/iJ 4 +SME 

Distance (pc) 396 ± 20 YY+a/i?*+SME 



a SME — "Spectroscopy Made Easy" package for the analysis of 
high-resolution spectra (Valcnti & Piskunov 1996). These parame- 
ters rely primarily on SME, but have a small dependence also on 
the iterative analysis incorporating the isochronc search and global 
modeling of the data, as described in the text. 

b YY+a/_R„+SME = Based on the YY isochrones flYi et aLII200lT) . 
a/R* as a luminosity indicator, and the SME results. 

Following iTorres et al.l (|2007f ). we explored the possi- 
bility that the measured radial velocities are not real, 
but are instead caused by distortions in the spectral line 
profiles due to c ontamination from a nearby unresolved 
eclipsing binary (jQueloz et al.l l2001). A bisector analy- 
sis b ased on the Keck s pectra was done as described in 
§5 of iBakos et al.1 (j2007) . We detect no correlated varia- 
tion between the bisector spans and the radial velocities 
(see Figure [5]). All but two of the bisector measurements 
are consistent with no variation. Fo llowi ng the meth- 
ods de scribed in lHartman et al.l (|2009T > and lKovacs et al.l 
(2010) we estimated the expected effect of contamina- 
tion from scattered moonlight on the bisectors, finding 
that the two outlier measurements correspond to the two 
spectra that are expected to be the most affected by 
sky contamination, when the Moon is 70° from the tar- 
get. Therefore, we conclude that the velocity variations 
are real, and that the star is orbited by a close-in giant 
planet. 

3.3. Global modeling of the data 

This section briefly describes the procedure we fol- 
lowed to model the HATNet photometry, the follow- 
up photometry, and the radial velocities simultaneously. 
More detail s on t he fitting methods can be found in 
IBakos et al.l ([2010). Our model for the follow-up light 
curves used the analytic formulae of iMandel fc Agoll 
(2002) with quadratic limb darken i ng coe fficients inter- 
polated from the tables bv iClaretl (|2004f ). The transit 
shape was parametrized by the ratio-of-radii p = i? p /i?*, 
the square of the impact parameter b 2 , and the recip- 
rocal of the half duration of the transit (/R*. We de- 



note this fitting set as {p,b 2 ,(/R*}- This set is cho- 
sen because of their simple geometric mean ings and the 
fact t hat these exhibit low correlations (see IBakos et al.l 
12010(1 . Our model for the H ATNet data was the si mpli- 
fied "P1P3" version of the IMandel fc Agol (12001 an- 
alytic functions following the method of IBakos et al.1 
( 2010). Following the formalism presented bv lPall (12009). 
the RVs were fitted with an eccentric Keplerian model 
parametrized by the semi- amplitude K and Lagrangian 
elements k = ecoscj and h = esincj, in which w is the 
longitude of periastron. 

Assuming a linear ephemeris, we assign the transit 
number N tr = to the complete follow-up light curve 
gathered on 2010 Feb 15. The eight main parameters 
describing the physical model were thus T c ^250, 7c,0: 
R p /R±, b 2 , (/-R*, K, k = ecosw, and h = esinw. 
Five additional parameters were included that have to 
do with the instrumental configuration. These are the 
HATNet blend factors i?i ns t,3i5 and i?inst,3i6, which ac- 
count for possible dilution of the transit in the HAT- 
Net light curves from background stars due to the broad 
PSF (20" FWHM), the HATNet out-of-transit magni- 
tudes -Mo, HATNet, 315 and M ,HATNot,3i6, and the relative 
zero-point 7 rc i of the Keck RVs. 

We extended our physical model with an instrumen- 
tal model that describes brightness variations caused by 
systematic e r rors i n the measurements as described in 
IBakos et all (|2010h . The HATNet photometry has al- 
ready been EPD- and TFA-corrected before the global 
modeling, so we only considered corrections for system- 
atics in the follow-up light curves. We chose the "ELTG" 
method, i.e., EPD was performed in "local" mode with 
EPD coefficients defined for each night, and TEA was 
performed in "global" mode using the sam e set of stars 
and T FA coefficients for all nights, as done in Bak os et al.l 

(poToh . 

The joint fit was accomplished us ing downhill sim- 
plex (AMOEBA;see iPress et~ail 119921 ) a nd the Ma rkov 
Chain Monte-Carlo method (MC MC. see lFord||2006l) us- 
ing "Hyperplane-CLLS" chains (B akos et al.l l2010h and 
the a naly t ic par tial derivatives for the transit light curve 
fr om I Pall (120091). A detailed description can be found 
in IBakos et al.l (|2010l) . The resulting geometric parame- 
ters pertaining to the light curves and velocity curves are 
listed in Table [5j Quotes values are the median and the 
error on the median from the a posteriori distribution of 
each parameter. 

Included in this table is the RV "jitter", which we 
added in quadrature to the internal errors for the RVs 
in order to achieve x 2 /dof = 1 from the RV data for 
the global fit. Auxiliary parameters not listed in the 
table are: T c _ 250 = 2454405.00857 ± 0.00162 (BJD), 
T c o = 2455243.81859 ±0.00029 (BJD), the blending fac- 
tors B in str,3i4 = 0.83 ± 0.05 and 5^,315 = 0.69 ± 0.03, 
and 7 rc i = 3.20 ± 2.62 ms _1 . The latter quantity rep- 
resents an arbitrary offset for the Keck RVs, and does 
not correspond to the true center of mass velocity of the 
system, which was listed earlier in Table [4] (7rv)- 

The planetary parameters and their uncertainties can 
be derived by combining the a posteriori distributions for 
the stellar, light curve, and RV parameters. In this way 
we find a mass for the planet of M p = 0.685 ± 0.033 Mj 
and a radius of R p = 1.242 ± 0.067 leading to a 
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mean density p p = 0.44 ± 0.07 gcm~ 3 . These and other 
planetary parameters are listed at the bottom of Ta- 
ble^ We note that the system may be slightly eccentric: 
e = 0.067 ± 0.024, u = 197 ± 36°. 

4. COMPARISON OF FITTING METHODS 

4.1. Alternative fitting parameter sets 

The fitting method adopted for the values quoted 
in Table [5] uses the parameter set {p£/R*,b 2 }, cho- 
sen for their low inter-parameter correlations. (,/R* 
is the reciprocal of the half-duration as computed us- 
ing an a pproximate express i on fo r the durati on com- 
ing fr om iTinglev fc Sackettl (|2005f ). Recently, iKippind 
(2010) showed that an improved approximate formula is 
possible for the duration. Replacing (,/R* with the recip- 
rocal of the new expression for the half-duration therefore 
offers gr e ater a ccuracy for the duration determination. 
IKippind (2010) labeled this parameter as T/i?» to be 
distinct from (,/R*. The new parameter set yields the 
greatest improvements for near-grazing, low-eccentricity 
orbits. 

To investigate the robustness of the results against 
different parameter sets, we refitted the EPD, TFA 
corrected HAT and FLWO light curves in conjunction 
with the radial velocities using the {p 2 ,6 2 ,T/i?»} pa- 
rameter set. Radial velocities are fitted using the La- 
grang e parameters as in Section 13.31 A new YY (|Yi et al.1 
120011 ) isochrone analysis is performed as described in 
Section 13.11 to show the effect on the physical parame- 
ters. The results are shown in the first results column of 
Table © 

The 6 2 parameter also gives us some pause for thought. 
One natural reason to select b 2 is that the duration of a 
transit is completely described in terms of b 2 i.e. there 
is no case of an isolated b which is not squared in the 
expression for the duration. Therefore, b 2 seems to be a 
natural parameter of the transit light curve. Further, b 2 
usually shows lower inter-parameter correlations than b. 

However, there are also two reasons why one should not 
choose b 2 . Firstly, geometrically b is more likely to have 
a uniform prior than b 2 . Therefore, low signal-to-noise 
transits will be biased towards higher impact parameters 
by fitting for b 2 . Secondly, b 2 cannot be negative and 
thus the posterior distribution of b 2 tends to get offset to 
an artificially more positive value due to the boundary 
condition that b 2 > 0. One resolution to this is to use 
b and let the parameter explore both negative and pos- 
itive values. Whilst a negative impact parameter may 
seem unphysical, the issue of its physicality is also irrel- 
evant since the transit is completely described in terms 
of b 2 and thus a negative b is always multiplied by it- 
self when computing the light curve morphology. Using 
b in this way permits for a symmetric distribution about 
6 = and thus improved estimates of the associated un- 
certainty. 

To investigate the value of fitting for 6, we repeated our 
global fits using the parameter set {p 2 , 6,21.5,3. 5}. Ti. 5,3.5 
replaces T/R* since although it may be slightly more 
correlated, it is more reasonable to expect a uniform prior 
on the duration than a uniform prior on its reciprocal. 
The results of these fits are shown in the third results 
column of Table |H 

The system values quoted between {p 2 ,b 2 , 



{p 2 ,b 2 , T/i?»} and {j> 2 ,6,2i.5,3. 5 } show excellent agree- 
ment and thus indicate that the system parameters are 
insensitive to the choice of priors. 

4.2. Fitted limb darkening 

We repeated the fits for the {p 2 ,b 2 , Y/i?*} and 
{p 2 , 6,21.5,3. 5} parameter sets with fitted limb darkening 
(second and fourth result columns of Table [6] respec- 
tively). In order to achieve convergence, we choose to 
only use linear limb darkening and thus fix the quadratic 
coefficient to be zero. 

The results show slight differences with the fixed limb 
darkening analogues. The most noticeable effect is in- 
creased error bars. Fitting for limb darkening means that 
our results are no longer dependent upon a stellar atmo- 
sphere model prediction, which is a clear desideratum. 
Our preferred final values are given by the {p ,6,21.5,3.5} 
parameter set with fitted limb darkening (last column of 
Table [51 with highlighted header). 

Both free limb darkening fits converge at u\ = 0.25 ± 
0.04. For 0.8 < \i < 1, the limb darkening from both 
the theoretical quadratic coefficients and the fitted coef- 
ficients is approximately equivalent. However, for low /1, 
the theoretical coefficients predict much stronger dark- 
ening effects than observed. Fitting a linear law through 
the quadratic coefficients gives u\ ln = 0.46, demonstrat- 
ing the stronger limb darkening predicted from theory. 

5. ORBITAL ECCENTRICITY 
5.1. Significance of the eccentric fit 

The eccentric fit suggests an orbital eccentricity of 
e = 0.052j:°;° 22 . In addition to the four fits dis- 
played in Table |6j we repeated our preferred model fit 
(i.e. {p 2 , 6,21.5,3.5} with fitted limb darkening) for a cir- 
cular orbit. This was done to provide a x 2 value for both 
fits, which can be used to infer the statistical significance 
of the eccentric fit. 

To accomplish this, we only take the x 2 from the radial 
velocity data, which dominates the determination of the 
Lagrangian orbital parameters and so the number of data 
points is n = 18. The period and time of transit are 
dominated by the photometry and very weakly affected 
by the few RV points. Therefore, these two degrees of 
freedom can be considered fixed. This leaves us with four 
degrees of freedom for an eccentric fit (k, h, K, 7) and 
two for the circular fit (K, 7). 

In evaluating the significance of the eccentric fit over 
the circular model it is important to penalize the ec- 
centric fit for using two extra degrees of freedom. We 
therefore choose to perform an F-test between the two 
models. The circular orbit fit has \ 2 = 59.5 and the 
eccentric fit has x 2 = 39.6. The false alarm probability 
from an F-test is eva luated to be 5.8% or 1.9 -q. We also 
performed the test of lLucv fc Sweeney! (|197lD . where the 
statistical significance of the eccentric fit is given by: 



P(e > 0) = 1 - exp 



e 

2^! 



(1) 



Where e is the modal value of the eccentricity, which 
is well approximated by the median for a uni mode dis- 
tribution. Using the ILucv &: Sweenevl (|1971l ) test, we 
find an eccentric fit is accepted at the 2.6-<r level. The 
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TABLE 5 
Orbital and planetary parameters 



Parameter Value 

Light curve parameters 

P (days) 3.355240 ± 0.000007 

T c (BJD) a 2455216.97667 ± 0.00028 

T M (days) a 0.1539 ± 0.0008 

Ti2 = T 34 (days) a 0.0141 ± 0.0006 

a/R* 7.58 ±0.35 

C/-R* 14.31 ± 0.06 

Rp/R* 0.0970 ± 0.0012 

>r o.o3«i •;;;!; 'v 

b = acosi/R* 0.189tog|o 

i (deg) 88.6 ± 0.7 

Limb-darkening coefficients b 

(Li (linear term) 0.1858 

bi (quadratic term) 0.3625 

RV parameters 

K (ms" 1 ) 83.0 ±3.4 

fc RV c -0.053 ±0.021 

/i RV c -0.017 ±0.042 

e 0.067 ±0.024 

co (deg) 197 ± 36 

RV jitter (ms -1 ) 7.4 

Secondary eclipse parameters 

T s (BJD) 2455218.542 ± 0.045 

T s 14 0.1492 ±0.0121 

T s ,i2 0.0137 ±0.0013 

Planetary parameters 

M p (Mj) 0.685 ±0.033 

R p (Rj) 1.242 ±0.067 

C(M p ,Rp) d 0.31 

p p (gem" 3 ) 0.44 ±0.07 

\ogg p (cgs) 3.04 ±0.05 

a (AU) 0.0465 ± 0.0006 

T cq (K) 1637 ± 42 

C 0.043 ±0.003 

Fper (lOSergs^cm- 2 ) 1.84 ±0.208 

F ap (^ergs-^m" 2 ) 1.44 ±0.168 

(F) (lOSergs" 1 cm" 2 ) f 1.62 ±0.169 

a T c : Reference epoch of mid transit that minimizes the correla- 
tion with the orbital period. It corresponds to Nt r — —8. BJD 
is calculated from UTC. T14: total transit duration. T12 — T34: 
ingress/egress time. 

Values for a quadratic law, adopted from the tabulations by 
Claret (2004) according to the spectroscopic (SME) parameters 
listed in Table g] 

c Lagrangian orbital parameters derived from the global modeling, 
and primarily determined by the RV data. 

Correlation coefficient between the planetary mass M v and ra- 
dius R p . 

e The Safronov number is given by — ^ ( V^sc/Vorb ) 2 — 
(a/R p )(M p /M„) (see|H arisen & Barman 2007). 

Incoming flux per unit surface area, averaged over the orbit. 

slightly higher significance likely comes from the fact this Tidal dissipation causes planetary orbits to circularize 
test does not penalize an eccentric model for using more over time. The maximum eccentricity a planet could 
degrees of freedom, whereas the F-test does. Therefore, initially have is e ~ 1. After N = 1 circularization 
based upon the current data for this system, an eccentric timescales, denoted r C i rc , the planet's eccentricity will re- 
fit is probable but not conclusive. duce by one e-fold, i.e. a factor of 2.72. For the planet to 
A resolution would be to obtain a secondary eclipse now have an eccentricity of e, the number of circulariza- 
measurement for the system, for which the mid-eclipse tion timescales which have transpired must be < — log(e) 
time would be dependent upon e cos w. The A: component and therefore TA gc < — log(e)r c i rc . This therefore con- 
dominates the eccentricity budget and thus its determi- strains the circularization timescale to be: 
nation would strongly constrain the eccentricity of this 

system. . > ^A ge ,„s 

clrc ^-log(e) [> 

5.2. Circularization timescales Using the method of | Adamg fc La J ((2QQ6I) , the 
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TABLE 6 
Orbital and planetary parameters 



Parameter 



{b 2 ,T/R*}, fixed LD {b 2 ,T/R*}, fitted LD {6,21.5,3.5}, fixed LD {b,Ti. 5 ,3. 5 }, fitted LD 



Light curve parameters 

^(days) 3.3552459t™°p 

T c (BJDutc - 2,450,000) a . . . . 5216.976721 ,; , , , 25 , 

Ti, 4 (s) a 13307.31^? ° 

Ti.5,3.5 (s) a 12112.1^4:o 

T 2 , 3 (s) a 10917.0l^| 

Ti, 2 ~ T 3 , 4 (s) a 1184.4l^'2 

(R P /R*) 2 (%) 0.93631™ 

II, ■ II- 0.09676lg-°°°« 

a/R* 7.7ll°;| 2 

T/B, (days" 1 ) 14.267±g:gg 

b n.oo •!;;;!;: 

0.000 •»;!>;';, 

i(deg) 90.01°;° 

/<• (gem" 3 ) 0.7701S 

Limb-darkening coefficients b 

ui (linear term) 0.1858* 

U2 (quadratic term) 0.3625* 

RV derived parameters 

fc R .v c -0.037l°'°" 

^Rv c -0-OlslS 

* 1 - UL " -0.004 

K (ms- 1 ) 82.6I32 

e 0.052 

w (deg) 206lg 

RV jitter (ms" 1 ) 7.53l 7 5°j 

log 9p (cgs) 3.0551°,;°^ 

Secondary eclipse parameters 

T s (BJDutc - 2,450,000) 5218.5751°,- °l° 

2Xi,4 (s) 128501^°° 
Stellar parameters 

M* (M ) 1-1841°;^ 

R* (Re) i.293±° °^ 

log(g*) (cgs) 4.287l°;°it 

L* (L Q ) 2.471^9 

M v (mag) 3.781°-^ 

Age (Gyr) 2.75l°;g 

Distance (pc) 404 _23 

Planetary parameters 

Mp (Mj) 0.680l°-°|° 

Rp (Rj) 1-2171™ 

C{Mp/R P } d 0.264 

Pr (gcm-3) 0.467™ 

» (AU) 0.046391™ 2 

T cq (K) 1623l|| 

(~> c 0761 + (ul0:i2 

w U - U ' D± -0.0()31 

F pcI (lOSergs-icm- 2 ) f l-™t°oil 

F ap (lO^rgs-icm- 2 ) f 1.43±g^| 

(F) (10 9 ergs- 1 cm -2 ) f 1-571^1 
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+ 3.53 
7.43 



5218.575l° ; ° 3 ° 
127301^° 



1 1«R+ 042 
1.18b_ .o4i 

i-294iS:SS 

4.286^°;°- 
2.48l°; 32 

q -Jj + 0.14 

,5 - 78 -0.14 
„ --+0.58 
*• ' a -0.65 
27 



405 



n fisi +° 031 

U.O»l_o.030 
1 O4Q + 072 

0.290 
439+0 069 

u.ioa 0.O66 

04641+° 00054 

U.U4O11_ 00054 

16241 44 
0745+° 0035 

U.Uf^O 0.0037 
.+0.20 
1 -0.14 
j + 0.16 
-0.18 



1.74] 
1.43^ 



1.57 



+0.18 
-0.15 



a T c : Reference epoch of mid transit that minimizes the correlation with the orbital period. It corresponds to Nt r — —8. T XjV : transit duration 
between contact points x and y. (x,y) — (1.5, 3.5) correspond to the sky-projected center of the planet overlapping the stellar limb, 
k Values for a quadratic law and fixed coefficients, adopted from the tabulations by Claret (2004) according to the spectroscopic (SME) 
parameters listed in Table [4] 

c Lagrangian orbital parameters derived from the global modeling, and primarily determined by the RV data. ^: Reciprocal of the 

scaling factor by which the true stellar density is modified from that found assuming a circular orbit (Kipping 2010). 

Correlation coefficient between the planetary mass M v and radius _R P . 
G The Safronov number is given by = 5 ( Vosc/ V orb ) 2 = {a/ R P )(M P / M t ) fsee lHansen fe Barmanll2007r) . 
Incoming flux per unit surface area, averaged over the orbit. 
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circularization timescale may be expressed as a function 
of the planet's tidal dissipative constant, Qp, for a low 
eccentricity system. 

r circ = Qp--— (--j (1-e) (3) 

Taking e ~ in the above expression together with 
the posterior distributions of the various parameters 
given above allows us to constrain Qp to be Qp > 
(6-1 jj o) x 10 6 . Note that the value above becomes 
eve n larger if we include th e (1 — e 2 ) 13 / 2 term from 
the I Adams &: Laughlinl (|2006t ) equation, but this requires 
some assumption of the histor y of the system. Given 
than Jupiter has Q P ~ 30,000 (jLainev et aT]|2009D . this 
limit raises some questions about why HAT-P-24b has 
such a lar ge value, somewhat similar to the situation for 
GJ 436b (jDeming et alJl2007[ ). In none of our 10 5 real- 
izations do we have a Qp value below 150,000 and there- 
fore in the absence of any eccentricity pumping, a large 
Q p value is a possible origin for the non-zero eccentricity. 
This is consistent with the ob servation of large Q p va lues 
in many other known TEPs (Matsumur a et alJ l2008). 

We note that for planets with initial eccentricities 
> 0.2, the above approximate expressions will be in- 
valid and a full backwards integration of the planet's 
orbital evolut i on wi ll be necessary, as pointed out by 
iLeconte et al.l (|2010| ). Such a detailed analysis remains 
outside of the scope of this paper, but our calculations 
do flag this system as possibly retaining an anomalously 
large eccentricity requiring further investigation. 

6. LINEAR RV DRIFT 

The unfolded residuals of an eccentric fit seem to hint 
at a negative linear drift in the radial velocities. We 
re-executed the global fit of the data including a RV gra- 
dient term j. We choose to use the {p 2 ,&,Ti. 5,3.5} param- 
eter set with fitted limb darkening again. The fits obtain 
7 = — 0.040^o;Q 2 g m/s/day with a slightly decreased ec- 
centricity of e = 0.048±^ 22 . By the ILucv fc Sweenevl 
(|197lD test, the eccentricity is now significant at the 2.3- 
a confidence level. We also note that the 7 parameter 
appears to have converge d in the MCMC trials w ith the 
Gelman- Rubins statistic (jGelman fc Rubinl Il992h satis- 
fying the criteria of being < 1 (value was 0.59), indicative 
of good-mixing. 

6.1. Statistical significance 
6.1.1. F-test 

There are numerous tests which one can employ to 
evaluate the significance of the gradient. The first one 
we tried was to compute the F-test between the eccentric 
orbit and the eccentric orbit + linear drift model. Pe- 
nalizing for one extra degree of freedom, the F-test find 
the drift model is accepted with 76.4% confidence. 

6.1.2. Odds ratio 

The second test we tried was to extract the posterior 
distribution of the gradient from the MCMC runs. If 
the gradient was equal to zero, we would expect 50% 
of the MCMC runs to give a positive value and 50% to 
give a negative value. In the eccentric + drift model, 
/ nog = 92.7% of the MCMC runs gave a negative 7. The 



odds ratio of the negative valued model over the 50:50 
model is: 

n _ °- 5 (a\ 

^drift/static — -, _ r ( 4 J 

L Jncg 

For only two possible models (i.e. a drift or static), the 
probability of the drift model being the correct one is 
P(drift) = 1 - [1/(1 + O dri ft/static)] = 87.2%. Therefore, 
both tests so far indicate a ~ 20% false alarm probability 
for the drift model. 

6.1.3. Bayesian Information Criterion 

The final test we performed was to re-fit all of the 
data using four possible models, using the {p 2 , b, Ti.5,3.5} 
and fitted limb darkening method, each with a different 
number of degrees of freedom, d: 

1. Circular orbit (d = 2) 

2. Circular orbit + linear drift (d = 3) 

3. Eccentric orbit (d = 4) 

4. Eccentric orbit + linear drift (d = 5) 

In each case we compute the Bayesia n Infor mation Cri- 
terion (BIC) (jSchwarzl (fl978h : ILiddld pOO^h ). given by 
BIC= x 2 + dlogn where n is the number of RV data 
points. BIC severely penalizes models for having more 
parameters and offers a statistically valid tool for model 
selection. We also compute the reduced x 2 , given by 
Xreduced = X 2 / ( n ~ c 0- As an example, the circular orbit 
has only two degrees of freedom in 7 and K. One might 
argue that P and tc are also degrees of freedom but in 
a global fit, which includes the HAT and FLWO time se- 
ries, these two parameters are overwhelmingly driven by 
the photometry and not the RV and thus the RV actually 
has negligible freedom in these parameters. 

The BIC model selection test indicates that the ec- 
centric orbit + linear drift model is the accepted model 
description of the current radial velocities when globally 
fitted with the current photometry for this system. We 
note that i) the eccentric models are consistently pre- 
ferred over the circular orbit models ii) the eccentricity 
is affected by a negligible degree by including the drift. 

6.1.4. Conclusion 

The linear drift model is the preferred model using the 
Bayesian Information Criterion. Other tests indicate the 
model is accepted with a false alarm probability of 20%. 
This is not sufficient to yet claim the trend is real and 
thus we encourage observers to obtain more observations 
to confirm or reject the existence of this trend. 

6.2. Properties of Putative HAT-P-24c 

We proceed here to constrain the properties of HAT-P- 
24c under the assumption the trend is real. Whilst this 
may not turn out to be true, it is useful to consider what 
the properties of the outer planet would be should the 
trend be later confirmed. 

The period of the outer planet would have to be much 
greater than the timescale of the observations, or a sinu- 
soidal pattern would have emerged and so P c ^> 265 days, 
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TABLE 7 

Comparison of four different models f or the RVs of 
hat-p-24b, described in section f6. 1 .31 



Model 


d 


•^reduced 


BIC 


7 


e 










m s~ 1 day — 




Circular 


2 


3.72 


65.3 


Q a 


a 


Circular + Drift . 


3 


3.35 


59.0 


-n n5^+ 026 

u.uoo_ 026 


ga 


Eccentric 


1 


2.83 


51.2 




o n52+ 022 


Eccentric + Drift 


5 


2.72 


49.9 


u.u*u_ 028 


048+ 022 



a Parameter is fixed. 

most likely of O ~ 1000 days which constrains a c > 2 AU 
by Kepler's Third Law. We note that the habitable 
zone pushes out as ahab = \J L*/Lq AU and occurs at 
~ 1.6 AU for HAT-P-24 and so the outer p lanet would 
likely be "cold". Using equation (1) from I Winn et al.l 
(2009), the gradient corresponds to an outer planet sat- 
isfying: 

= M = 0.0821°;^ M,AU- 2 (5) 

a c (_r 

Based upon the a c constraint, this therefore implies 
M c sini c > 0.3 Mj, Aside from the RVs, there are 
observational consequences for HAT-P-24b due to the 
outer planet. The system will behave as an inner- 
outer binary and thus the outer planet will induce a 
light-time travel effect, potentially detectable with tran- 
sit timing variations (TTV) of HAT-P-24b. The inner 
binary will orbit the barycentre with semi-major axis 
{a c M c sin i c )/(M* + Mb). Therefore the peak-to-peak 
light-time effect, for an outer planet of negligible eccen- 
tricity, will be: 

TTVught (peak - to - peak) ~ Ui -% (6) 

Adding in the best-fit value of 7 gives 
TTV ligh t(peak-to-pcak) = 0.078(P c /years) 2 s. 
Therefore, we require P c > 3.6 years for a > 1 second 
TTV and would need P c > 36 years for a > 100 s TTV. 
We also evaluated the TTV effect due to a distant, 
pertur bing planet, as described in case IV of lAgol et al.l 
pOOl . We find that a 1000 day period 0.3 Mj would 
generate an r.m.s. TTV of 0.008 s for e c = 0.1, 0.06 s 
for e c = 0.5 and 0.8 s for e c = 0.9. The challenge 
of measuring TTVs to this precision over such long 
time-scales is a daunting one and unlikely to reap any 
reward with current instrumentation. 

7. TRANSIT TIMES 

Whilst the photometric quality of the FLWO data is 
sufficient for full free fitting, three of the four FLWO 
light curves are only partial transits and so the errors on 
the duration and therefore mid-transit time diverge for 
unconstrained fitting parameters. 

A solution to this is to work under the assumption 
that the duration and depth of the transit do not change 
from transit-to-transit. As a result of the partial transits, 
a transit duration variation (TDV) analysis is therefore 
not possible but TTVs can be obtained provided it is 
understood that they are derived as inherently model 



TABLE 8 
Fitted mid-transit times. 



Epoch 


tc/(BJD UT c-2,450,000) 
days 


O-C 
seconds 


-11 
-9 
-3 



5206.910851™ 2 
5213.62078t™ 5 5 

koto 7K01 7 +0.00080 
OZOO.(OOl/_ Q0082 

5243.81878+°" 


-fi q+ 53 - 6 
u ' y -53.6 

-55 fi+ 56 - 2 

90 1+69.1 
zo - ± -70.8 
11 Q+26.0 
1J - a _26.0 



dependent values, where the model is that of constant 
duration and depth. 

We therefore extract the parameters p 2 , b, T, esinw, 
ecosw, P and the linear limb darkening coefficient u\ 
from the posterior distribution of the global fit. We se- 
lect the {b,T} parameter set with fitted limb darkening 
as our favorite solution for this purpose. The free param- 
eters of the individual transit fits are OOT and tc- We 
stress that the parameters assumed to be constant from 
transit-to-transit are still allowed to float around their 
median- value with standard deviation given by their de- 
rived uncertainties. This ensures the errors are correctly 
propagated into the mid-times. The final times are given 
in Table M 

Using the linear ephemeris derived from the global fit, 
including all HAT data, we find that the FLWO transits 
show no excess variance yielding \ 2 — 1.3 for 2 degrees of 
freedom. The RMS of the four O-C values is 34.7 seconds. 
With only four transit times, it is not possible to conduct 
a meaningful TTV analysis. However, these transits may 
be used a benchmark for future TTV searches on this 
system. 

8. FOLLOW-UP POSSIBILITIES 

8.1. Secondary Eclipse 

The Keck RVs indicate that e cos lu is significantly non- 
zero and thus suggests the secondary eclipse of HAT-P- 
24b would occur with a timing offset from that a c i rcular 
orbit. Neglecting terms of order cot 2 i, iSternd ([1940D 
shows that the timing offset is given by: 

PI e cosw-\/l — e 2 / ecosw \ \ 

At = —\ -— -r^r + arctan - (7 

7T y(l -(e sin uj) 2 ) \VT^Jj 

Using the posterior distribution of parameters from the 
{p 2 ,6,Ti. 5i3 . 5 } fitted-limb-darkening MCMC run, we es- 
timate that the secondary eclipse should occur (1.9 ± 
0.7) hours earlier than that expected for a purely circu- 
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lar orbit. The detection of the eclipse would therefore 
strongly constrain e cos u> which would allow for a revised 
global fit to the data. Constraining ecosw in this way 
would allow us to re-compute the statistical significance 
of the RV linear drift model over the static model. 

For Spitzer's 3.6 /im and 4.5 /zm channels, assuming 
uniform redistribution of energy around the planet and 
zero albedo, we estimate a depth of ~ 0.075% and 
~ 0.10%. At K s = 10.5 mag, the brightness is compara- 
ble to that of TrES-4 (K s = 10.3 mag) which has been 
observed with Spitzer and eclipses measured to preci- 
sions o f 0.011% and . 016% for the two channels respec- 
tively (|Knutson et al J 120091) . What HAT-P-24b lacks in 
a slighter fainter host star it makes up for with a slightly 
lon ger duration than TrES- 4b (~ 12000 s for TrES-4b, 
see iMandushev et all ([2057h . and 13150 s for HAT-P- 
24b). We therefore estimate that the secondary eclipses 
will be detectable with SNRs of ~ 7 and ~ 6 for 3.6^m 
and 4.5/im respectively. 

8.2. Transmission Spectroscopy 

Molecular constituents in the terminator of the atmo- 
sphere can absorb light and cause the transit depth to 
increase. The spectral variations in the transit depth al- 
low for the detection of molecules within an exoplanet's 
atmosphere. This depth change can be estimated from 
first principles by calculating the scale height of the at- 
mosphere and computing the e xpected chan ge in transit 
depth using equation (36) from|Wjnn| (|2010j ). 

Using the posterior distribution of parameters from 
the {b,Ti 535} fitted-limb-darkening MCMC run, we es- 
timate that AS ~ (0.014 ± 0.001)iV ff % where N H is 
the number of scale heights of atmosphere absorbed by 
the molecular species (of order unity) , and we have used 
jiM = 2a.m.u for H^. Using the same Spitzer uncer- 
tainty estimates from before, we would require Nh ~ 3 
for even H 2 to be detectable. Therefore, HAT-P-24b 
would likely be a challenging target for transmission 
spectroscopy. 

8.3. Rossiter-McLaughlin Effect 

Owing to HAT-P-24's relatively ra pid stellar rota - 
tion, a large Rossit er-McLaughlin effect (Ros siterl <| 1924f t : 
McLaughlin (1924)) amplitude is expected. Equally, a 
large sig nal for spectral line t omography could be de- 
tect able_JCameron^FaL| I2010D . Using equation (40) 
from iWinnl (|2010f ). and the posterior distribution of 
parameters from the {6,11.5,3.5} fitted-limb-darkening 
MCMC run, we estimate AV RM = (95 ± 5) m/s. Given 
the RV measurements are essentially jitter-limited at 
7.2 m/s, this indicates that we expect a very large 
signal-to-noise for the RM effect of HAT-P-24b, reach- 
ing SNR~13. 

Another motivation for measuring the planet's RM ef- 
fect is that HAT-P-24 has an effective temperature of 
T off = (6373 ± 80) K and thus lies above the ~ 6250 K 
threshold for wh ich most system s eem to exhibit signifi- 
cant obliquities (|Winn et al.ll2010D . This therefore indi- 
cates that we can not only expect a very large RM effect 



b ut also possibly a high ly oblique configuration. 

iGaudi fc Winfil ([2007) showed that v sin i and A be- 
come degenerate for low impact parameter transits and 
so the large RM amplitude predicted for this system 
should help in solving for the system parameters. 
9. SUMMARY 

We announce the detection of a 0.68 Mj transiting 
exoplanet on a 3.36 day orbit around an F8 star (sys- 
tem parameters are found in the last column of Ta- 
ble O . We find that the planet retains a small eccen- 
tricity of e = 0.052l ; 017 with a 5.8% false alarm prob- 
ability, which may suggest either a perturbing planet in 
the system or low tidal dissipation within the planet of 
Qp > (6 ± 3) x 10 6 . Most of the eccentricity originates 
from the ecosu term and thus we predict that a sec- 
ondary eclipse observation, which is shown to be feasible, 
should confirm/reject the eccentricity hypothesis conclu- 
sively. 

We have performed a detailed investigation of the ef- 
fects upon the system parameters by using different fit- 
ting sets. In three different parameter sets, we find a 
consistent solution indicating the result is not sensitive 
to the priors. The effects of fixing versus fitting limb 
darkening coefficients are also investigated, which leads 
to slightly increased error bars but a consistent best-fit 
solution. 

Using the Bayesian Information Criterion as a model 
selection tool, we find the Keck radial velocities are 
best described by a model consisting of non-zero or- 
bital eccentricity and a negative linear drift of (—14.6 ± 
10.2) m/s/year with a false alarm probability of 20%. We 
consider this trend to be currently not statistically sig- 
nificant, but warranting further investigation. 

HAT-P-24 has a relatively rapid stellar rotation of 
i>sinz = (10.0 ± 0.5)km/s, and we therefore predict 
HAT-P-24b should exhibit one of the largest known 
Rossiter-McLaughlin effect amplitudes for an exoplanet 
(AVrm — 95 m/s) . F u rther it has recently been sug- 
gested by I Winn et"aT1 (|2010l ) that hot-stars have com- 
panions on preferentially oblique orbits and so HAT-P- 
24 (T cff = (6373 ± 80) K) would be an excellent target to 
further investigate this hypothesis. 
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